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I. INTRODUCTION 



Modern tactical missiles require increased stand-off ranges 
and need to meet threats from highly maneuverable air targets. 

The high maneuverability of air targets has directed the use of 
defence missiles capable to develop higher lift accelerations and 
more complex control laws. In order to accomplish the 
requirement of large stand-off ranges, propulsion systems using 
air-breathing engines have been studied and developed in recent 
years. The advent of air-breathing engines has naturally led to 
the consideration of BTT missiles in order to minimize the inlet 
angle-of-attack . 

The necessity of more complex control laws has introduced the 
application of modern control and estimation theory, since more 
complicated information of the missiles states are needed. BTT 
controlled missiles are generally characterized by increased 
maneuverability and considerable drag reduction over conventional 
cruciform, roll stabilized STT controlled missiles. Certain 
limitations in technology [Ref. 1] have delayed the development 
of BTT control systems, and consequently any progress in the area 
of BTT autopilots. 

Major technological improvements during the last decade, as 
the availability of advanced digital computers, reopened the 
issue and made BTT control feasible in spite the added complexity 
of control laws for the autopilots. In addition, certain types 
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of ramjet engines [Ref. 2], which are candidate propulsion 
systems for modern tactical mission requirements of range and 
high altitude [Ref. 3], have presented a need for a missile control 
technique to maintain effective inlet flow. This was the main 
reason for given further impetus to the investigation and 
development of BTT control. 

Despite the fact that BTT steering may provide improved 
performance for a missile system, there are still unanswered 
questions concerning stability during homing phase, guidance 
performance, autopilot guidance logic and subsystem requirements. 
All these questions have to be investigated and properly answered 
in order for BTT steering to be considered as a viable control 
method for high performance missiles. 

During the past decade many missile programs [Ref. 3] were 
initiated to improve the capability of steering tactical missile 
via BTT control with results that have greatly advanced the 
understanding of the various missile subsystems. In the 
autopilot area many different types have been designed and 
developed. All of them force the missile to roll or bank, so that 
the steering maneuver occurs with the missile axis oriented in a 
specific or preferred direction with respect to the incoming 
airstream. This class of autopilots is usually known as POC 
autopilots . 

The main criterion for the selection of a particular type of 
autopilot is based upon the guidance, airframe and propulsion 
system requirements. Generally, missiles with either one or two 
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planes of symmetry use a POC autopilot which forces the missile 
to bank in order to turn as an aircraft. If this motion is 
coordinated, then the autopilot is referred to as a CBTT 
autopilot . 

In the guidance area, radome aberation effects for frequency 
guidance are of major concern [Ref. 3] and are being investigated 
in great extent. Also, the interaction between BTT control, 
antenna stabilization and sensor orientation are some of the 
additional concerns that have to be properly addressed. However, 
simplified studies [Ref. 4], which neglect radome effects and 
assume that the missile motion is entirely coordinated, have 
proven that BTT control can provide acceptable performance with 
roll rates that are not excessive for autopilot design. These 
studies were made for a medium range area and long suppression 
mission, and considered both high lift (i.e., planar) and 
moderate lift (i.e., cruciform) airframe configurations. 

In order to take full advantage of CBTT control, planar 
airframes have been designed to increase the lifting capability 
in -one direction without the weight and drag penalty associated 
with orthogonal lifting surfaces [Ref. 5]. These airframes have 
aerodynamic properties characterized by increased potential to 
enhance CBTT control. 

The present thesis addresses the design and analysis of 
discrete lateral autopilots for application to CBTT missiles. 

The first part reviewed the design procedure of the two 
individual lateral channels as developed by Arrow [Ref. 6]. The 
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design was performed using classical techniques and involved the 
uncoupled yaw and roll channels for the elliptical and circular 
airframes respectively. The resulted continuous open loop 
designs were analyzed in terms of their transient and frequency 
responses and found to be in accordance with the desired 
requirements specified in [Ref. 6], Then, applying analog-to- 
digital conversion, the corresponding discrete lateral autopilots 
were obtained and analyzed. 

The second part utilized modern control design methods to the 
already discussed discrete single- input lateral autopilots. This 
allowed comparison with the preceding classical design, and more 
importantly established a technique to extend some of the results 
to the more general multivariable case. At first, assuming 
availability of all states for feedback purposes, application of 
the Ackermann formula led to a discrete state-feedback designed 
autopilot. Next, since the state vector of the state-feedback 
model is not usually accessible to direct measurement, an estimator 
was introduced as an additional dynamic design in order to 
implement control to the original system. The state-feedback and 
estimator designed autopilots were analyzed for both lateral 
channels and found to have satisfactory responses. 

Finally, coupling the discrete pitch and roll channel 
autopilots, the state-feedback and estimator designs were 
obtained and proved to be robust. 

The analysis in all the above cases was performed using the 
existed at Naval Postgraduate School OPTSYS and ORACLS Fortran 
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program for the continuous and discrete systems respectively. 
Additionally, the last part that dealt with the coupled pitch and 
roll channel autopilot utilized the POPLAR design program 
developed by Gordon [Ref. 7]. 
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II. CLASSICAL DESIGN AND ANALYSIS OF LINEAR UNCOUPLED 

LATERAL AUTOPILOTS 



A. GENERAL 

The initial phase in the design of lateral CBTT autopilots 
involved the design and analysis of the individual uncoupled 
lateral channels (i.e., yaw and roll) with prescribed 
relationships between speeds-of-response. These relationships 
when coupled with the corresponding ones of the longitudinal 
uncoupled channel (i.e., pitch) would meet the requirements of 
the overall CBTT autopilot. 

The uncoupled autopilot design method was classical and used 
a combination of frequency response and root locus techniques. 
Utilization of this particular design method led to the 
achievement of practical bandwidths (i.e., sufficient high 
frequency attenuation), and. in turn provided the range of required 
missile body angular rates and control motions. In addition, the 
resulting design minimized the influence of aerodynamic 
variations on desired' responses. The application of the 
uncoupled channels to the whole CBTT autopilot was accomplished 
by an appropriate choice of the relative time constants of the 
individual channels. In order to achieve the desired maneuver 
plane acceleration the roll channel was designed to have a time 
constant of 0.5 seconds. The yaw uncoupled channel, which 
follows the roll motion to produce the required coordination 
(i.e., minimization of sideslip angle), is designed to have a 
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more rapid response with a time constant of 0.39 seconds for the 
circular airframe. The detailed requirements for the classical 
design of the uncoupled yaw and roll channel autopilots are 
presented in Appendix A. 

A fixed flight condition (i.e., constant Mach number and 
altitude) was selected for this preliminary performance study. 
Fixed flight conditions are typically used in autopilot designs 
to identify and cure critical areas of concern. When autopilot 
requirements are satisfied at fixed flight conditions, then areas 
of concern caused by varying them are addressed. The selected 
flight condition of 60000 feet altitude and Mach number 3.95 
provided sufficient dynamic pressure, so that the missile 
maneuvers resulted in large enough angles-of-attack to exercise 
sideslip control. Aerodynamic data for this particular flight 
condition are provided in Appendix A. 

The aerodynamic models developed for stability studies in the 
frequency domain were linearized about a trim angle-of-attack for 
both lateral channels. The following three assumptions were 
made : 



1. The plane Xg - -E& of Figure 2.1 was the maneuver plane. 

2. The missile was trimmed in pitch (i.e., M y = 0, at fixed 
values of -X , q, and 

Rather than use the assumption that the missile roll rate is 
approximately zero as it is normally done for the roll stabilized 
STT control, the following assumption was made for BTT: 
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Figure 2.1 Aerodynamic Sign Convention and Axis 



3. The missile roll rate was constant. 

Linearized aerodynamic derivatives are given in Appendix C. 

In this chapter the analysis of both continuous and discrete 
uncoupled lateral channels of a BTT autopilot was based on the 
transient and frequency responses of maneuver plane accelerations, 
body angular rates and tail incidence angles. A general block 
diagram of a BTT autopilot with all its channels is shown in 
Figure 2.2. Inertial acceleration command were applied in polar 
coordinates (i.e., magnitude of the command ft applied to pitch 
and the direction to roll autopilot). The yaw autopilot was 

slaved to the roll autopilot in order to minimize the yaw and 
roll motions. Achieved maneuver plane accelerations in 
rectangular coordinates (i.e., and fty ) were determined by 

resolving achieved body-fixed accelerations (i.e., J] ^ and f) y . ) 
through missile roll rate (i.e., Euler angles 0 and ijJ were 
assumed to be sufficiently small). 

B. AIRFRAME CONFIGURATIONS 

The two airframe configurations studied in this work were 
taken from [Ref. 6] and are shown in Figure 2.3 and Figure 2.4. 
Although the configuration in Figure 2.3 reveals a body of 
circular cross section and that of Figure 2.4 an elliptical one, 
both airframes have the same cross sectional area distribution. 

In specific, the circular cross sectional body has a closure 
ratio Abase/Amax of 0.69 with Amax occurring at 68% missile body, 
whereas the elliptical airframe has a 3:1 cross section. 
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Figure 2.2 Bank-to-Turn Autopilot 
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Figure 2.3 Model of Circular Airframe Configuration I /6-scale 
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Both airframe configurations are tail-controlled using four 
identical control surfaces which are located flush with the body- 
base with a +30° dihedral. In the case of the elliptical body, 
the hinge line was skewed such that a 10° control deflection 
measured at the body-tail juncture had a resultant 7.04° surface 
deflection. Thus, the aerodynamic control effectiveness in terms 
of deflection measured at the body-tail is lower for the 
elliptical airframe although it is nearly the same in terms of 
resultant surface deflection. 

The total span of the mono-wings is the same for each 
configuration, which results in larger wing area for the circular 
airframe. The wing area and span for the circular airframe were 
chosen as typical of current maneuvering missiles. The wing for 
the elliptical concept was determined by projecting the 
elliptical body on the circular body-wing planform. The resultant 
when exposed wing planform became the wing for the elliptical 
body . 

Comparison of the elliptical airframe with the corresponding 
circular indicates the following: 

1. About 30$ more normal force that is nearly independent of 
angle-of-attack can be achieved at supersonic speeds. 

2. Values of longitudinal stability Darameter C_ are more 

‘I 

positive, and with more pronounced nonlinearities in 
pitching moment at subsonic speeds. 

3. Levels of directional stability are increased and more 
compatible with levels of longitudinal stability. 
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4. More yaw control is available although suitable locations 
for tails on the body are more limited because of the 
geometry of the elliptical airframe. 

The two airframe configurations were sized to provide 
realistic geometric and mass properties. The details are 
presented in Appendix D. 

C. UNCOUPLED YAW CHANNEL AUTOPILOT FOR ELLIPTICAL AIRFRAME 

The purpose of the uncoupled yaw channel autopilot of a CBTT 
missile is to minimize the sideslip angle ( j3 ) , or provide 
coordinate motion between the yaw and roll channels. The easiest 
way to accomplish this is by designing the uncoupled yaw channel 
(i.e., roll and pitch dynamic effects neglected) as a regulator 
(i.e., no guidance command and with rate and acceleration 
feedback) to help minimize the sideslip angle. 

A block diagram of the uncoupled yaw channel is shown in 
Figure 2.5. In this diagram both the aerodynamic model and yaw 
control law are involved. The normal acceleration ( }] v ) is not 

I c 

used to command the CBTT autopilot. Instead, it is used for the 
design and analysis of the uncoupled channel. The command used by 
the coupled system is shown in dashed lines and is a yaw angular 
rate command (r Q ). The yaw control law shown in Figure 2.6 [Ref. 
6] is governed by missile body angular rate (r) and yaw normal 
acceleration (j]y). At the flight condition of interest (i.e., 

60 kft altitude and Mach number 3.95) the yaw control law 
determines the required command ( «r v ) to an actuator which is 
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Figure 2.5 Uncoupled Yaw Channel 
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Figure 2.6 Yaw Control Law 



approximated as a first-order lag at 30 Hz. The rate compensator 
computes the high frequency attenuation and is used to minimize 
aerodynamic variations on the quality of the regulator. On the 
other hand, the acceleration compensator measures the acceleration 
bandwidth via the time constant of the acceleration response of • 

I / 

1 . Transfer Functions of Aerodynamic Model 

The aerodynamic transfer functions of the uncoupled yaw 
channel autopilot are: 



^ t %£) fcC-AE + 6 c ) ^ * * 
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(II .C. 1-6 ) 



- _ (57.5) 

t zz 

/B54- 

N r y (II. C. 1-7) 



Substituting the values of aerodynamic data 
linearized aerodynamic derivatives (Table VIII), and 
and mass properties (Table IX), equations (II. C. 1-3) 
( II . C . 1 -7 ) become : 

A = -0.1351 
5 = 0.0436 
C = -17.2748 
E = 34.5495 
K = 0.4823 



(Table VII), 

geometric 

through 



(II. C. 1-8) 
(II. C. 1-9) 
(II .C. 1-10) 
(II. C. 1-1 1 ) 
(II .C. 1-12) 



Introducing the above equations (II. C. 1-8) through 
(II. C. 1-12) to (II. C. 1-1) and (II. C. 1-2), the aerodynamic 
transfer functions of the uncoupled yaw channel for the circular 
airframe and zero angle-of-at tack can be obtained. 

a. Transfer Function of Yaw Angular Rate: 
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(II. C. 1-13) 



b. Transfer Function of Yaw Normal Acceleration: 
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(II. C. 1-14) 
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Rearranging the above transfer functions in terms of the 



variable pairs (r, Cy ) and (n , , ) respectively, they become: 



zsi 0. 0444 Z$ - 1S..J8I 5 r » - 56. €6 S T s -X. 74-2 2S y 


(II .C. 1-15) 
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(II. C. 1-16) 



Applying inverse Laplace transformation, the following 
set of linear differential equations can be obtained: 



z iO'04-94- z -IS. 2S15 c=-56>66 by - 2. 7422 by 
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Both equations (II. C. 1-17) and (II.C.1-18) form a second 

order system of linear differential equations in which the 

forcing function involv'es derivative terms. Using rules of 

state-space representation [Ref. 8] the following equations are 

obtained : 

» 
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2 . 



Equations of Yaw Control Law and Actuator 



a. Acceleration Compensator Equation 

The acceleration compensator equation obtained from 



Figure 2.6 is: 

= 0. 31 <746 
y ” CU 5tl 




(II. C. 2-1 ) 



Rearranging and applying inverse Laplace transforma- 
tion, (II. C. 2-1) turns into the following linear differential 
equation : 

y - ~ 5 y n y - 1 . 5975 n Y ^ (n.c. 2 - 2 ) 

Substituting equation (II. C. 1-25) into (II. C. 2-2), 
the last becomes: 

y - V 5975 2^ - 3. ?534-3y - 1.5975 ily^ (II .c.2-3 ) 

b. Rate Compensator Equation 

The rate compensator equation also obtained from 



Figure 2.6 is: 

- _ 4.85 (4r SH ) 



c'ytz) 



(II. C. 2-4) 



Rearranging and applying inverse Laplace transforma- 
tion to (II. C. 2-5), it turns into: 

o = 0.485^ + 0.485 r f 4. S5y+- 4-. S5 jt (II. C. 2-5) 

Substituting equations (II. C. 1-23), (II. C. 1-24) and 



(II. C.2-3) into the above, it becomes: 
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Actuator Equation 

The actuator equation obtained from Figure 2.5 is: 



c . 




(II. C. 2-7) 



Rearranging and applying inverse Laplace transforma- 



tion, the last equation turns into: 



r ' 1 98 . 4 - c - 1 98 , 4 Oy 



(II. C. 2-8) 



3 . Design Approach and Analysis of Continuous System 
Utilizing state-space representation, the equations 
(II. C. 1-19) through (II. C. 1-22), (II. C. 2-2) (II. C. 2-6) and 
(II. C. 2-8) can be modeled in a seventh-order system of the form 
x=Fx+Gu. The continuous plant system and input matrixes F and G 
are shown in Table I, and the state vector is: 




(II .C. 3-1 ) 



where the state variables are: 



x-| , x£ : yaw angular rates 

z-j, Z 2 *. yaw normal accelerations 



y 



output of acceleration compensator 



input command in the actuator 



yaw tail incidence 
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PLANT SYSTEM AND INPUT MATRICES; UNCOUPLED YAW 
CHANNEL AUTOPILOT; CLASSICAL DESIGN; CONTINUOUS OPEN 
LOOP SYSTEM; ELLIPTICAL AIRERAME 
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Executing the OPTSYS program, using an input step 
function which represents "1 gee command" at zero trim angle-of 
attack, the pole-zero and time and frequency response plots are 



obtained . 

The pole-zero plot 
continuous open loop system is 
are : 

S ; - ~ 1 74-. 3- 16 

5ji =■ - 6. \ 0^2. -h j 1C . 6 596 

5 $ ~~ 6.1C4-2 -j 10.6596 

S 4 = -0.02 /747s 4-.0S919 
5 S - -G.C2 17475- j 4CS919 
S 6 : -2.9296 
S 7 : -2.9296 - j 199929 



of Figure 2.7 indicates that the 
table, since the s-plane poles 

(II .C. 3-2) 
(II. C. 3-3) 
(II. C. 3-4) 
(II. C. 3-5) 
(II. C. 3-6) 
(II. C. 3-7) 
(II .C. 3-8) 



The time response plots of the yaw normal acceleration, angular 
rate and tail incidence are shown in Figures 2.8 through 2.10. 

In particular the yaw normal acceleration time response plot has 
a 0.39 seconds time constant, 1% overshoot and a steady-state 
error of 0.018. These results are in accordance with the 
requirements referred in Appendix A, that is a time constant of 
0.4 seconds, overshoot less than 10% and steady-state error not 
necessarily equal to zero. All the above three time response 
plots are identical with those presented in [Ref. 6]. 

Figures 2.11 through 2.16 show the frequency response 
plots of the yaw normal acceleration, angular rate and tail 
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incidence, from which the phase crossover frequencies and gain 
margins of Table II can be obtained. The positive gain and phase 
margins of the open loop system ensure the relative stability of 
the closed loop (controlled) system. 

4 . Design Approach and Analysis of Discrete System 

Utilizing analog- to-digital conversion by the aid of 
ORACLS program and for a sample period of 0.0125 seconds, a 
seventh-order discrete system of the form x(k+1 ) = Ax(k)+Bu(k) is 
obtained. The discrete pant system and input matrices A. and B 
are shown in Table III. 

The pole-zero plot of Figure 2.17 indicates that the 
discrete open loop system is also stable, since the z-plane 



poles are: 

Zj = C. II5G/6' (II.C.4-1) 

C. 916554- tj O, \228£Z (II.C.4-2) 

Zx - C.9i g'554-j 0, 122.86 Z . (II.C.4-3) 

Z ^ •= 6. 996 422 F j O.Ob \Q (II. C. 4-4) 

Z S * -j O.OS \C7SS ( II . C . 4-5 ) 

5 90S SC-5 +j 0,0 (II.C.4-6) 

Z 7 =0.<?05S6s-jO.O561546 , tt r , 7 , 



The time response plots of the yaw normal acceleration, 
angular rate and tail incidence for the discrete uncoupled yaw 
channel are presented in Figures 2.18 through 2.20. A close 
observation of these plots indicates that they are identical with 
those of the continuous classical system found in the previous 
section . 
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Figure 2.7 Pole-Zero Plot; Uncoupled Yaw Channel Autopilot; 

Classical Design Continuous Open Loop System; 
Elliptical Airframe 
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YAW NORMAL ACCELERATION, HY 
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Figure 2.8 Yaw Normal Acceleration vs Time; Uncoupled Yaw 
Channel Autopilot; Classical Design; Continuous 
Open Loop System; Elliptical Airframe 
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YAW ANGULAR RATE.R 
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Figure 2.9 Yaw Angular Rate vs Time; Uncoupled Yaw Channel 
Autopilot; Classical Design; Continuous Open 
Loop System; Elliptical Airframe 
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YAW TAIL INCIDENCE, <5Y 
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Figure 2.10 Yaw Tail Incidence vs Time; Uncoupled Yaw Channel 

Autopilot; Classical Design; Continuous Open 
Loop System; Elliptical Airframe 
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Figure 2.11 Yaw Normal Acceleration-Gain vs Frequency; Uncoupled 

Yaw Channel Autopilot; Classical Design; Continuous 
Open Loop System; Elliptical Airframe 
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PHASE - DEGREES 
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Figure 2.12 Yaw Normal Acceleration-Phase vs Frequency; 

Uncoupled Yaw Channel Autopilot; Classical 
Design; Continuous Open Loop System; Elliptical 
Air f rame 
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MAGNITUDE - DECIBELS 
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Figure 2.13 Ya w Angular Rate-Gain vs Frequency; Uncoupled 

Yaw Channel Autopilot; Classical Design; 
Continuous Open Loop System; Elliptical 
Airframe 
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Figure 2.14 Yaw Angular Rate-Phase vs Frequency; Uncoupled 

Yaw Channel Autopilot; Classical Design; 
Continuous Open Loop System; Elliptical 
Airframe 
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MAGNITUDE - DECIBELS 
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Figure 2.15 Yaw Tail Incidence-Gain vs Frequency; Uncoupled 

Yaw Channel Autopilot; Classical Design; 
Continuous Open Loop System; Elliptical 
Airframe 



51 



PHASE - DEGREES 

116.0 - 94.0 - 72.0 - 50.0 - 28.0 - 6.0 16.0 38.0 60.0 82.0 




Figure 2.16 Yaw Tail Incidence-Phase vs Frequency; Uncoupled 

Yaw Channel Autopilot; Classical Design; 
Continuous Open Loop System; Elliptical 
Air f rame 
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PHASE CROSSOVER FREQUENCIES AND GAIN MARGINS; 
UNCOUPLED YAW CHANNEL AUTOPILOT; CLASSICAL DESIGN; 
CONTINUOUS OPEN LOOP SYSTEM; ELLIPTICAL AIRFRAME 
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D. UNCOUPLED ROLL CHANNEL AUTOPILOT FOR CIRCULAR AIRFRAME 

The uncoupled roll channel autopilot of a CBTT missile is 
commanded to roll the missile so as to put the preferred maneuver 
direction in the direction of the guidance acceleration command. 
The desired maneuver plane acceleration should be attained as 
rapidly as the achieved body-fixed pitch acceleration. To 
accomplish this, the uncoupled roll channel autopilot (i.e., yaw 
and roll dynamic effects neglected) was designed to have the roll 
angle time constant equal to the time constant of the normal 
acceleration achieved by the uncoupled pitch channel autopilot. 

A block diagram of the uncoupled roll channel is shown in 
Figure 2.21. In this diagram both the aerodynamic model and roll 
control law are involved. The roll control law shown in Figure 
2.22 is commanded by roll angle (<p c ) and governed by roll angular 
rate (p) and roll angle ((jl). 

The design and analysis of the uncoupled roll channel 
autopilot was performed in this section for the stable at zero 
angle-of-attack circular airframe. 

1 . Transfer Functions of Aerodynamic Model 

The aerodynamic transfer functions of the uncoupled roll 
channel autopilot obtained from Figure 2.21 are: 

a. Transfer Function of Roll Angular Rate: 

1 A Cj «- 

^ _ x v < C' £ r- 

I ~ j (II. D. 1-1) 

° x XX 
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PLANT SYSTEM AND INPUT MATRICES; UNCOUPLED YAIV 
CHANNEL AUTOPILOT; CLASSICAL DESIGN; DISCRETE OPEN 
LOOP SYSTEM; ELLIPTICAL AIRPRAME 
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Figure 2.17 Pole-Zero Plot; Uncoupled Yaw Channel Autopilot; 

Classical Design; Discrete Open Loop System; 
Elliptical Airframe 
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YAW NORMAL ACCELERATION, I1Y 

0.0 0.1 0.2 0.3 0.4 0.5 0.8 0.7 0.8 0.9 




Figure 2.18 Yaw Normal Acceleration vs Time; Uncoupled Yaw 

Channel Autopilot; Classical Design; Discrete 
Open Loop System; Elliptical Airframe 



57 



YAW ANGULAR RATE.R 

0.05 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.00 0.65 




Figure 2.19 Yaw Angular Rate vs Time; Uncoupled Yaw 

Channel Autopilot; Classical Design; 
Discrete Open Loop System; Elliptical 
Airframe 
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YAW TAIL INCIDENCE, <5Y 
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Figure 2.20 Yaw Tail Incidence vs Time; Uncoupled Yaw 

Channel Autopilot; Classical Design; 
Discrete Open Loop System; Elliptical 
Airframe 
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b. Transfer Function of Roll Angle: 

= _L p (II. D. 1-2) 

Substituting the values of aerodynamic data (Table VII), 
linearized aerodynamic derivatives (Table VIII) and geometric and 
mass properties (Table IX), equation (II. D. 1-1) becomes 

ps = 8.03462 (II. D. 1-3) 

Applying inverse Laplace transformation to equations 
(II. D. 1-2) and (II. D. 1-3), the following set of linear 
differential equations is obtained: 

p = 8.03462 o, k (II .D. 1-4) 

(p = p (II. D. 1-5) 

2 . Equations of Roll Control Law and Actuator 
a. Roll Angle Compensator Equation 

The roll angle compensator equation obtained from 



Figure 2.22 is: 

2. Z 

X - s * C <t c -q>) (II. D. 2-1) 

IT' 1 

Rearranging and applying inverse Laplace 
transformation (II. D. 2-1) turns into the following linear 
differential equation: 

X = -17. 6$ - 8X + 17.6 ^ (II. D. 2-2) 

b. Rate Compensator Equation 

The rate compensator equation also obtained from 



Figure 2.22 is: 




aisc«(— n) 



( x-p) 



(II .D. 2-3) 
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Figure 2.21 Uncoupled Roll Channel 
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Figure 2.22 RoJl Control Law 



Rearranging and applying inverse Laplace transforma- 
tion to (II. D. 2-3), it turns into: 

Y + 5Y = 0.0502X - 0.0502p + 0.75325X - 0.75325p (II. D. 2-4) 

Substituting equations (II. D. 1-4) and (II. D. 2-2) into 
(II. D. 2-4), the last becomes: 

Y=-0.7532p-0. 883526+0. 35 165X-5Y-0. 40334^ +0.88352% (II .D. 2-5) 
c. Pseudo-Differential Equation 

The oseudo-dif f er ent ia 1 eauation obtained from Fi sure 



2.22 is: 



*1* 



0. C! 5055 s 



P 



(II .D .2-6) 



Rearranging and applying inverse Laplace 
transformation to (II. D. 2-6), it turns into: 

X 1 = -6X 1 + 0.0781 98p (II. D. 2-7) 

Substituting equation (II. D. 1-4) into (II. D. 2-7), the 
last becomes: 

X 1 = - 6 X 1 + 0.628291 <5^ (II. D. 2-8) 

d. Equation of Actuator Compensator 

The equation of actuator compensator obtained from 
Figure 2.22 is: 

S r fi 

■ C Y-Xj) (II.D.2-9) 



°^ s 



iO 



jr + i 
1 b 



tion 




Rearranging and applying inverse 

to (II.D.2-9), it turns into: 

» • 

=-15 % c +0.13636Y-0.13636X 1 +15Y-15X 1 



Laplace 



transf orma- 



(II .D. 2-10) 
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Substituting equations (II. D. 2-5) and (II. D. 2-8) into 
(II. D. 2-10), the last becomes: 

--CUCV7/ p -aUC45tp-l-6).047' r /5X>i4:5iS2r-l4:(S/fX 1 - Cj2L4Sq ! 

(III .D. 2-1 1 ) 



e. Actuator Equation 

The actuator equation obtained from Figure 2.21 is: 




57.5 




(II .D. 2-12) 



Rearranging and applying inverse Laplace 
transformation, the last equation turns into: 



^ = 10795.32 -188.4 



(II. D. 2-13) 



3 . Design Approach and Analysis of Continuous System 
Utilizing state-space representation, the equations 



(II. D. 1-4), (II. D. 1-5), (II. D. 2-2), (II. D. 2-5), (II. D. 2-8), 
(II.D.2-11), and (II. D. 2-13) can be modeled in a seventh-order 
system of the form x=Fx+Gu. The continuous plant system and 
input matrices F and G are shown in Table IV, and the state 
vector is: 



XCK) = 



P 

X 

r 



(II. D. 3-1 ) 
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where the state variables are: 



p : roll angular rate 

(j) : roll angle 

X : output of roll angle compensator 
Y : output or rate compensator 

X : output of pseudo-differentiator 

6^ : input command in the actuator 

: roll tail incidence 

Executing the OPTSYS program, using an input step 
function which represents "1 gee command" at zero trim angle-of- 
attack, the pole-zero and time and frequency response plots are 
obtained . 

The pole-zero plot of Figure 2.23 indicates that the 
continuous open loop system is stable, since the s-plane poles 
are : 

5 = - 174.785 (roll angular rate) (II. D. 3-2) 

Sj = -9. 25097+ j 28. 4098 (roll angle) (II. D. 3-3) 

% = -9 . 25097- j 2 8 . 4098 (output of roll angle compensator) 

(II. D. 3-4) 

- -2 . 46608+ j2 . 7 1 1 52 (output of rate compensator network) 

(II. D. 3-5) 

S ^ = -2 . 46608- j2 . 7 1 1 52 (output of pseudo-differentiator) 

(II. D. 3-6) 

Sg - -8.98209 (input command in the actuator) (II. D. 3-7) 

Sj = -5.19932 (roll tail incidence) (II. D. 3-8) 

The time response plots of the roll angle, angular rate 
and tail incidence are shown in Figure 2.24 through 2.26. In 
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TABLE IV 
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particular, the roll angle time response plot has a 0.55 seconds 
time constant, 3 % overshoot and a steady-state equal to zero. 
These results are in accordance with the requirements referred in 
Appendix A, that is a time constant of 0.5 seconds, overshoot 
less than 10% and zero steady-state roll angle error. All the 
above three time response plots are identical with those 
presented in [Ref. 6]. 

Figures 2.27 through 2.32 show the frequency response 
plots of the roll angle, angular rate and tail incidence, from 
which the phase crossover frequencies and gain margins of Table 
V can be obtained. The positive gain and phase margins of the 
open loop system ensure the relative stability of the closed loop 
(controlled) system. 
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Figure 2.23 Pole-Zero Plot; Uncoupled Roll Channel Autopilot 

Classical Design; Continuous Open Loop System; 
Circular Airframe 
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ROLL ANGLE , <? 

0.0 0.1 0.2 0.3 0.4 0.5 0.« 0.7 0.8 0.9 




Figure 2.2k Roll Angle vs Time; Uncoupled Roll Channel 

Autopilot; Classical Design; Continuous Open 
Loop System; Circular Airframe 



69 




ROLL ANGULAR RATE.P 

0.0 O.I 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 




Figure 2.25 Roll Angular Rate vs Time; Uncoupled Roll Channel 

Autopilot; Classical Design; Continuous Open Loop 
System; Circular Airframe 
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ROLL TAIL INCIDENCE, <5R 
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Figure 2.26 Roll Tail Incidence vs Time; Uncoupled Roll Channel 

Autopilot; Classical Design; Continuous Open Loop 
System; Circular Airframe 
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MAGNITUDE - DECIBELS 
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Figure 2.27 Roll Angle-Gain vs Frequency; Uncoupled Roll Channel 

Autopilot; Classical Design; Continuous Open Loop 
System; Circular Airframe 
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PHASE - DEGREES 
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Figure 2.28 Roll Angle-Phase vs Frequency; Uncoupled Roll 

Channel Autopilot; Classical Design; Continuous 
Open Loop System; Circular Airframe 
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MAGNITUDE - DECIBELS 

80.0 - 62.0 -44 0 -26 0 - 8.0 10.0 26.0 46.0 64.0 02.0 




Figure 2.29 Roll Angular Rate-Gain vs Frequency; Uncoupled Roll 

Channel Autopilot; Classical Design; Continuous Open 
Loop System Circular Airframe 
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PHASE - DEGREES 
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Figure 2.30 Roll Angular Rate-Phase vs Frequency; Uncoupled Roll 

Channel Autopilot; Classical Design; Continuous Open 
Loop System; Circular Airframe 
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MAGNITUDE - DECIBELS 
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Figure 2.31 Roll Tail Incidence-Gain vs Frequency;. Uncoupled 

Roll Channel Autopilot; Classical Design; 
Continuous Open Loop System; Circular Airframe 



76 



PHASE - DEGREES 

336.0 - 300.0 - 264.0 - 228.0 - 102.0 - 156.0 - 120.0 - 84.0 - 48.0 - 12.0 







Figure 2.32 Roll Tail Incidence-Phase vs Frequency; Uncoupled 

Roll Channel Autopilot; Classical Design; 
Continuous Open Loop System; Circular Airframe 
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TABLE V 



PHASE CROSSOVER FREQUENCIES AND GAIN MARGINS; UNCOUPLED ROLL 
CHANNEL AUTOPILOT; CLASSICAL DESIGN CONTINUOUS OPEN LOOP SYSTEM; 

CIRCULAR AIRFRAME 







PHASE CROSSOVER 
FREQUENCY (rad/sec) 


GAIN 

MARGIN (db) 


ROLL 


ANGLE (d>) 


10.0461 


60.2165 


ROLL 


ANGULAR RATE (p) 


10.0461 


60.2165 


ROLL 


TAIL INCIDENCE (^ ) 


— 


— 
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4 . Design Approach and Analysis of Discrete System 



Utilizing analog- to-digital conversion by the aid of 
ORACLS program and for a sample period of 0.0125 seconds, a 
seventh-order discrete system of the form x(K+1 ) = Ax(k)+ Bu(k) is 
obtained. The discrete plant system and input matrices A and B 
are shown in Table VI. 

The pole-zero plot of Figure 2.33 indicates that the 
discrete open loop system is also stable, since the z-plane poles 
are : 

- 0.0992805 (roll angular rate) (II. D. 4-1) 

z z - 0.835216+jO. 309736 (roll angle) (II. D. 4-2) 

- 0 . 8352 1 6 - jO . 309736 (output of roll angle compensator) 

(II .D. 4-3) 

■^ 4 - = 0 . 969087 +J 0 . 0328588 (output of rate compensator network) 

(II. D. 4-4) 

25 = 0 . 969087 - jO . 0328588 (output of pseudo-differentiator) 

(II. D. 4-5) 

2 ^. = 0.893797 (input command in the actuator) (II. D. 4-6) 

Zj - 0.937075 (roll tail incidence) (II. D. 4-7) 

The time response plots of the roll angle, angular rate 
and tail incidence for the discrete uncoupled roll channel are 
presented in Figures 2.3^ through 2 . 36 . A close observation of 
these plots indicates that they are identical with those of the 
continuous classical system found in the previous section. 
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Figure 2.33 Pole-Zero Plot; Uncoupled Roll Channel; Classical 

Design; Discrete Open Loop System; Circular Airframe 
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Figure 2.34 Roll Angle vs Time; Uncoupled Roll Channel 

Autopilot; Classical Design; Discrete Open Loop 
System ; 
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Figure 2.35 Roll Angular Rate vs Time; Uncoupled Roll Channel 

Autopilot; Classical Design; Discrete Open Loop 
System; Circular Airframe 
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Figure 2.36 Roll Tail Incidence vs Time; Uncoupled Roll Channel 

Autopilot; Classical Design; Discrete Open Loop 
System; Circular Airframe 
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III. MODERN CONTROL DESIGN AND ANALYSIS OF LINEAR 
UNCOUPLED LATERAL AUTOPILOTS 



A. GENERAL 

The task of this chapter is the design and analysis of the 
same discrete uncoupled lateral channel autopilots discussed in 
the previous chapter, using different techniques which are based 
on modern control formulation. The difference in the two 
approaches is entirely in the design method since the end result, 
a set of difference equations providing control, is identical. 

Modern control theory is contrasted with the classical 
control theory in that the former is applicable to multi-input- 
multi-output systems, which may be linear or nonlinear time- 
invariant or time-varying, while the latter is applicable only to 
linear time-invariant single-input-single-output systems. Also, 
modern control theory is essentially a time-domain approach, 
while the conventional classical control theory is a complex 
frequency-domain approach. 

System design in classical control theory is based on trial- 
and-error procedures which, in general, will not yield optimal 
control systems. System design in modern control theory, on the 
other hand, enables the design of optimal control systems of 
great complexity and good accuracy with respect to given 
performance indexes. In addition, design in modern control 
theory can be carried out for a class of inputs instead of 
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specific input function, such as the impulse, step or sinusoidal 
functions and can also include initial conditions. 

One of the most attractive features of modern control design 
method is that the procedure consists of two independent steps. 

One step assumes that all the system states are available for 
feedback purposes. In general, even if this is not a practical 
enough assumption since it needs a large number of sensors, it is 
usually adopted in order to accomplish the first design step, 
namely the control-law. The remaining step is the design of an 
estimator which estimates the entire state vector, given 
measurements of portion of the state provided by the system output 
equation. The final control algorithm consists of the control- 
law and estimator combined, where the control-law calculations 
are based on the estimated states rather than the actual states. 
This substitution is reasonable and the combined design can give 
closed loop characteristics which are unchanged from those 
assumed in designing the control-law and estimator separately. 



B. DISCRETE STATE-FEEDBACK DESIGN 



Considering the following discrete control system: 



X ( K+1 ) = Ax(k) + Bu ( k ) 
Y ( k ) = Hx ( k) 



(III .B. 1-1 ) 
(III .B. 1-2) 



the control-law design is also referred to as s ta te-f eedback 
design since it is simply the feedback of a linear combination of 
all the system states, that is: 



u ( k ) = Fx ( k ) 



(III .B. 1-3 ) 



where F: control-law gain vector 
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Thus, the characteristic equation of the controlled (closed loop) 
system is: 

det ( zI-A+BF ) = 0 (III. B. 1-4) 

The discrete state-f eedback design, providing that the system is 
controllable, consists then of finding the control-law gain 
vector F so that the roots of (III. B. 1-4) are in desirable 
locations . 

A program logic for computing the control-law gain vector via 
the Ackermann's formula is given in Appendix E [Ref. 9]. 

Utilizing this control algorithm a Fortran program was written 
(Appendix F) which has as inputs the sample period, the discrete 
plant system and control input matrices, the s-plane poles and 
provides as output the control-law gain vector. 

1 . Uncoupled Yaw Channel for Elliptical Airframe 
a. Control-Law Gain Vector 

Executing the Ackermann Fortran program described in 
Appendix F with inputs: 

(1) Sample period of 0.0125 seconds 

(2) Discrete plant system and control input matrices A and B 
of Table III 

(3) S-plane poles defined in equations (II. C. 3-2) through 
(II.C.3-8) 

the following control-law gain vector for the elliptical airframe 

of the uncoupled yaw channel is obtained: 

F = [ — 1 .01 95 -0.109 -0.3492 0.0393 -0.61 52 32.4775 -31.0868] 

(III .B. 1-5) 
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b. Design Approach and Analysis 

The discrete state- feedback designed yaw autopilot 
can be found by introducing the control-law gain vector of 
(II. B. 1-5) into the original system. 

The pole-zero plot of Figure 3.1 indicates that the 
discrete closed lop system is stable, since the z-plane closed 
loop poles are: 

2 ( = 0.46461 (yaw angular rate) (III. B. 1-6) 

= 0.898888+jO. 137039 (III. B. 1-7) 

2j = 0 . 898888- jO . 1 37039 (yaw normal acceleration) (III. B. 1-8) 

2 4 = 0.998423+jO. 051078 (III. B. 1-9) 

25 = 0.998423- jO. 05 1 0789 (output of acceleration compensator 

network) (III. B. 1-10) 

Z c - 0 . 962974+ jO . 03587 1 2 (input command in the actuator) 

(III.B. 1-11) 

zy = 0. 962974-j0. 0358712 (yaw tail incidence) (III.B. 1-12) 

The time response plots of the yaw normal 
acceleration, angular rate and tail incidence are presented in 
Figures 3.2 through 3.4. A close observation of the above pole- 
zero and time response plots for the discrete yaw state-f eedback 
design indicates that they are identical with those of the 
discrete classical design found in the previous chapter. 

c. Simplified Design 

The discrete state- feedback designed autopilot of the 
previous section can be simplified by reducing the returning gain 
loops. This can be accomplished by placing zeros into appropriate 
elements of the control-law gain vector: 
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F= [ - 1 .0 1 95 -0.109 -0.3^92 0 -0.6152 32.4775 -31.0868] 

(III. B. 1-13) 

The pole-zero and time response plots of the 
resulting simplified state-feedback yaw autopilot, shown in 
Figures 3.5 through 3.8, do not present significant differences 
from the corresponding plots of the discrete classical design 
apart from the overshoot which was slightly increased. The 
discrete closed loop system is again stable, since the z-plane 
closed loop poles are: 

Z, = 0.467984 ) III . B . 1-14) 

Z^ = 0.891328+j0. 126019 (III. B. 1-15) 

z 5 = 0.891328-jO. 126019 (III. B. 1-16) 

Z^. = 0.998422 + jO. 0510788 (III. B. 1-17) 

z 5 = 0.998422-jo. 0510788 (III. B. 1.18) 

Z c = 0.969075+jO. 0429695 (III. B. 1-19) 

Z- r = 0.969075-jO. 0429695 (III. B. 1-20) 

2 . Uncoupled Roll Channel for Circular Airframe 
a. Control-Law Gain Vector 

Following the same procedure as in the yaw channel 
case apart from the use of discrete matrices A and B from Table 
III and s-plane poles from (II. D. 3-2) through (II. D. 3-8), the 
control-law gain vector for the circular airframe roll channel 
was found to be: 

F=[0 0.0004 0 0.0178 -0.0179 -1.499 0.0023] (III. B. 2-1) 
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Figure 3.1 Pole-Zero Plot; Uncoupled Yaw Channel Autopilot; 

State-Feedback Design; Discrete Closed Loop System 
Elliptical Airframe 
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YAW NORMAL ACCELERATION, HY 
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Figure 3.2 Yaw Normal Acceleration vs Time; Uncoupled Yaw 

Channel Autopilot; State-Feedback Design; Discrete 
Closed Loop System; Elliptical Airframe 
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Figure 3*3 Yaw Angular Rate vs Time; Uncoupled Yaw Channel 
Autopilot; State-Feedback Design; Discrete 
Closed Loop System; Elliptical Airframe 
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YAW TAIL INCIDENCE,^' 
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Figure 3.4 Yaw Tail Incidence vs Time; Uncoupled Yaw Channel 
Autopilot; State-Feedback Design; Discrete 
Closed Loop System; Elliptical Airframe 
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Figure 3.5 Pole-Zero Plot; Uncoupled Yaw Channel Autopilot; 

Simplified State-Feedback Design; Discrete Closed 
Loop System; Elliptical Airframe 
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Figure 3.6 Yaw Normal Acceleration vs Time; Uncoupled Yaw 
Channel Autopilot; Simplified State-Feedback 
Design; Discrete Closed Loop System; Elliptical 
Airframe 
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Figure 3.7 Yaw Angular Rate vs Time; Uncoupled Yaw 
Channel Autopilot; Simplified State- 
Feedback Design; Discrete Closed Loop 
System; Elliptical Airframe 
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Figure 3-8 Yaw Tail Incidence vs Time; Uncoupled Yaw 
Channel Autopilot; Simplified State- 
Feedback Design; Discrete Closed Loop 
System; Elliptical Airframe 
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b. Design Approach and Analysis 

The discrete state- f eedback designed roll autopilot 
can be found by introducing the control-law gain vector of 
(III. B. 2-1) into the original system. 

The pole-zero plot of Figure 3*9 indicates that the 
discrete closed loop system is stable, since the z-plane closed 
loop poles are: 

2, = 0.0989714 (roll angular rate) (III. B. 2-2) 

Z Z = 0.832909+jO. 315088 (roll angle) (III. B. 2-3) 

= 0 . 832909- j0 . 31 5088 (output of roll angle compensator) 

(III. B. 2-4) 

2*. = 0 . 969236+ jO . 0333773 (output of rate compensator network) 

(III. B. 2-5) 

^5 = 0 . 969236 - jO . 0333773 (output of pseudo-differentiator) 

(III .B. 2-6 ) 

2^ = O .896583 (input command in the actuator) (III. B. 2-7) 

= 0.938233 (roll tail incidence) (III. B. 2-8) 

The time response plots of the roll angle, angular 
rate and tail incidence are presented in Figures 3.10 through 
3.12. A close observation of the above pole-zero and time 
response plots for the discrete roll state-feedback design 
indicates that they are identical with those of the discrete 
classical design found in the previous chapter. 

c. Simplified Design 

The state-feedback roll designed autopilot can be 
simplified by reducing the returning gain loops as follows: 

F=[0 0 0 0.178 -0.0179 -1.499 0] (III. B. 2-9) 



98 



The pole-zero and time response plots of the 
resulting simplified state- feedback roll autopilot, shown in 
Figures 3.13 through 3.16, do not present significant differences 
from the corresponding plots of the discrete classical design. 

The discrete closed loop system is again stable, since the z- 
plane closed loop poles are: 

Z, = 0.0992745 

= 0.836632+jO. 305719 
Z 5 = 0.836632-jO. 305719 
^4 = 0.967647+jO. 0334948 
% = 0.967647-jO. 0334948 
^6 = 0.895937 
2-7 = 0.938256 

C. DISCRETE ESTIMATOR DESIGN 

The state-feedback design discussed in the last section 
assumed that all system states were available for feedback 
purposes. Since the state vector is not always accessible to 
direct measurement, an estimator is going to be introduced in 
this section as an additional dynamic design in order to 
implement control to the original system. The estimator design 
method consists mainly of determining algorithms which will 
reconstruct all the states, given measurements of a portion of 
them. 



(III. B. 2-10) 
(III. B. 2-1 1 ) 
(III .B. 2-12) 
(III. B. 2-13) 
(III .B. 2-14) 
(III. B. 2-15) 
(III .B. 2-1 6 ) 
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Figure 3-9 Pole-Zero Plot; Uncoupled Roll Channel Autopilot; 

State-Feedback Design; Discrete Closed Loop System; 
Circular Airframe 
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Figure 3.10 Roll Angle vs Time; Uncoupled Roll Channel Autopilot; 

State-Feedback Design; Discrete Closed Loop System; 
Circular Airframe 
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Figure 3.11 Roll Angular Rate vs Time; Uncoupled Roll Channel 

Autopilot; State-Feedback Design; Discrete Closed 
Loop System; Circular Airframe 
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Figure 3.12 Roll Tail Incidence vs Time; Uncoupled Roll Channel 

Autopilot; State-Feedback Design; Discrete Closed 
Loop System; Circular Airframe 
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Figure 3.13 Pole-Zero Plot; Uncoupled Roll Channel Autopilot; 

Simplified State-Feedback Design; Discrete Closed 
Loop System; Circular Airframe 
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Figure 3.14 Roll Angle vs Time; Uncoupled Roll Channel Autopilot; 

Simplified State-Feedback Design; Discrete Closed 
Loop System; Circular Airframe 
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Figure 3.15 Roll Angular Rate vs Time; Uncoupled Roll Channel 

Autopilot; Simplified State-Feedback Design; 
Discrete Closed Loop System; Circular Airframe 
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Figure 3.16 Roll Tail Incidence vs Time; Uncoupled Roll Channel 

Autopilot; Simplified State-Feedback Design; 
Discrete Closed Loop System; Circular Airframe 
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Considering the same as in the state- f eedback case discrete 
control system, a prediction estimator defined by the following 
equation is introduced: 

x(K+1) = Ax (k) + Bu(k) + K[X(k) - Hx(k)] (III. C. 1-1) 

A 

where x: estimate state vector 
K: estimate gain vector 

In this closed loop estimator, shown in Figure 3.17, the 
difference between the measured and estimated output is fed back 
and the model is constantly corrected with this error signal 

A 

which is defined as x=x-x. The difference equation describing 
the behavior of the error is obtained by subtracting equation 
from the actual plant output equation (II. B. 1-2): 

x(K+1 ) = [A - KH] x(k) (III. C. 1-2) 

Thus the characteristic equation of the controlled (closed loop) 
system is: 

det (zl - A + KH) = 0 (III. C. 1-3) 

The discrete estimator design, providing that the system is 
observable, consists then of finding the estimator gain vector K 
so that the roots of (III. C. 1-3) are at desirable locations. 

The estimator gain vector can be obtained again, as in the 
case of the state-f eedback design, by application of the 
Ackermann program of Appendix F with inputs the sample period, 
the transposes of the discrete plant and output matrices, and 
faster s-plane poles from those of the continuous system. 
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Figure 3.17 Discrete Closed Loop Ustinia tor 



1 . Uncoupled Yaw Channel for Elliptical Airframe 



a. Estimator Gain Vector 

Executing the Ackermann Fortran program (Appendix F) 

with inputs: 

(1) Sample period of 0.0125 seconds 

(2) The transposes of the discrete yaw system plant and output 

T T 

matrices, that is A 1 and H-| , where: 




1 

0 

0 

0 

0 

0 

0 



(III.C. 1-4) 



(3) S-plane poles slightly 
open loop system, that 
5, = -174.386 

5^ = -6. 1 1 42+ j 1 0 . 6396 • 

S 6 - -6. 1 1 42- j 1 0 . 6396 
= -0.0227475+j4. 08919 
S 5 = -0.0227475-J4. 08919 
$6 = -2 . 9396+ j2 . 99929 
$T = -2 . 9396- j2 . 99929 
the transpose of the estimator 
airframe of the uncoupled yaw 
from which the following K is 



faster than those of the continuous 
is with more negative real parts: 

(III.C. 1-5) 
(III.C. 1-6) 
(III.C. 1-7) 
(III.C. 1-8) 
(III.C. 1-9) 
(III.C. 1-10) 
(III.C. 1-11) 

gain vector for the elliptical 
channel is calculated as output, 
obtained : 
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K = 



-0.0034 

0.0008 

-0.0003 

-0.0177 

0.0005 

-0.0017 

-0.0240 



(III .C. 1-12) 



b. Design Approach and Analysis 

The discrete estimator designed yaw autopilot can be 
found by introducing the estimator gain vector of (III. C. 1-12) 
into the original system. 

The pole-zero plot of Figure 3*18 indicates that the 
discrete closed loop system is stable, since the z-plane closed 
loop poles were found to be: 

2 j = 0.117115 (yaw angular rate) 

Z z = 0.918174+jO. 122926 

•£5 = 0 . 9 1 8 1 74- jO . 1 22926 (yaw normal acceleration) 

- 0.998399+jO. 0510958 



(III.C. 1-13) 
(III. C. 1-14) 
(III.C. 1-15) 
(III.C. 1-16) 



2 . - = 0 . 998399 - j 0. 051 0958 (output of acceleration compensator 



network) 



(III .C. 1 -17) 



- 0 . 963249+ jO . 036 1 306 (input command in the actuator) 

(III.C. 1-18) 

Z T = 0.963249-j0 . 0361306 (yaw tail incidence) (III.C. 1-19) 

The time response plots of the yaw normal 
acceleration, angular rate and tail incidence are presented in 
Figures 3.19 through 3-21. A close observation of the above pole- 
zero and time response plots of the discrete yaw estimator design 
indicates that they are very close to those of the discrete 
classical design found in the previous chapter. 
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Figure 3.18 Pole-Zero Plot; Uncoupled Yaw Channel Autopilot; 

Estimator Design; Discrete Closed Loop System; 
Elliptical Airframe 
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c. Simplified Design 

The discrete estimator designed autopilot of the 
previous section can be simplified by reducing the returning gain 
loops. This can be accomplished by placing zeros into 
appropriate elements of the estimator gain vector: 



K 



-0.0034 

0 

0 

-0.0177 

0 

-0.0017 

-0.0240 



(III .C. 1-20) 



The pole-zero and time response plots of the 
resulting simplified estimator yaw autopilot, shown in Figures 
3.22 through 3.25, do not present significant differences from 
the corresponding plots of the discrete classical design, apart 
from the overshoot which was slightly increased. The discrete 
closed loop system is again stable, since the z-plane closed loop 
poles are: 



2, =0.117117 

^ = 0.91 8434+ jO . 122898 
= 0. 91 8434-jO . 122898 
^4. = 0.998422+J0. 0510954 
^5 = 0.998422-j0. 0510954 
= 0.962965+j0. 0357874 
= 0. 962965-j0 .0357874 



(III .C. 1-21 ) 
(III. C. 1-22) 
(III .C. 1-23) 
(III.C. 1-24) 
(III .C. 1-25) 
(III.C. 1-26) 
(III.C. 1-27) 
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2 . 



Uncoupled Roll Channel for Circular Airframe 



a. Estimator Gain Vector 

Following the same procedure as in the yaw channel 
case apart from the use of: 

(1) The transpose of the discrete roll system plant matrix 

(2) S-Plane poles slightly faster than those of the continuous 
open loop roll system, that is: 



5, = -184.795 

5^ = -9.26097+J28.4098 

% = -9.26097-j28.4098 

S+ = -2 . 47608+ j2 .71152 
S 5 = -2.47608-j2.71 152 
^ = -8.99209 
S 7 = - 5.20032 



(III. C. 2-1 ) 
(III. C. 2-2) 
(III. C. 2-3) 
(III. C. 2-4) 
(III. C. 2-5) 
(III. C. 2-6) 
(III .C. 2-7) 



the estimator gain vector for the circular airframe roll channel 
is obtained: 



K 



f ound 
into 



0.004 

0 

0 

- 0.0002 

0.0003 

0 

-0.0492 
— — 

b. Design Approach and Analysis 

The discrete estimator designed roll 
by introducing the estimator gain vector 
the original system. 



(III. C. 2-8) 



autopilot can be 
of (III. C. 2-8) 
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Figure 3.19 Yaw Normal Acceleration vs Time; Uncoupled Yaw 

Channel Autopilot; Estimator Design; Discrete 
Closed Loop System; Elliptical Airframe 
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Figure 3.20 Yaw Angular Rate vs Time; Uncoupled Yaw 

Channel Autopilot; Estimator Design; 
Discrete Closed Loop System; Elliptical 
Air f rame 
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Figure 3.21 Yaw Tail Incidence vs Time; Uncoupled Yaw 

Channel Autopilot; Estimator Design; 
Discrete Closed Loop System; Elliptical 
Airframe 
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Figure 3.22 Pole-Zero Plot; Uncoupled Yaw Channel 

Autopilot; Simplified Estimator Design; 
Discrete Closed Loop System; Elliptical 
Airframe 
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Figure 3-23 Yaw Normal Acceleration vs Time; Uncoupled Yaw 

Channel Autopilot; Simplified Estimator Design; 
Discrete Closed Loop System; Elliptical Airframe 
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Figure 3.24 Yaw Angular Rate vs Time; Uncoupled Yaw 

Channel Autopilot; Simplified Estimator 
Design; Discrete Closed Loop System; 
Elliptical Airframe 
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Figure 3-25 Yaw Tail Incidence vs Time; Uncoupled Yaw 

Channel Autopilot; Simplified Estimator 
Design; Discrete Closed Loop System; 
Elliptical Airframe 
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The pole-zero plot of Figure 3.26 indicates that the 



discrete closed loop system is stable, since the z-plane closed 
loop poles were found to be: 





— 


0.0967149 


(III. C. 2-9) 


** 


- 


0.830876+jo. 31913 


(III. C. 2-10) 


Z 5 


r 


0.830876-jo. 31913 


(III .C. 2-1 1 ) 




- 


0.970259+jo. 0325175 


(III. c. 2-13) 


*5 


= 


0.970259-jO. 0325179 


(III. C. 2-14) 






0.894774 


(III. C. 2-15) 




= 


0.937225 


(III. C. 2-16) 






The time response plots of the 


roll angle, angular 



rate and tail incidence are presented in Figures 3*27 through 
3.29. A close observation of the above pole-zero and time 
response plots for the discrete roll estimator design indicates 
that they are very close to those of the discrete classical design 
found in the previous chapter. 

c. Simplified Design 

The estimator roll designed autopilot can be 
simplified by reducing the returning gain loops as follows: 



K 



0 

0 

0 

0.0178 
-0.0179 
-1 .499 
0 



(III .C. 2-17) 



The pole and time response plots of the resulting 
simplified estimator roll autopilot, shown in Figures 3-30 
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through 3-33, do not present significant differences from the 
corresponding plots of the discrete classical design. The 
discrete closed loop system is again stable, since the z-plane 
closed loop poles are: 



Z, = 0.0992745 

= 0.836632+jo. 305719 
= 0.836632-jO. 305719 
z 4 = 0.967647+j0. 0334948 
= 0.967647-jO. 0334948 
^ = 0.895937 

^7 = 0.938256 



(III .C. 2-1 8) 

(III. C. 2-19) 

(III .C. 2-20) 
(III. C. 2-21 ) 
(III .C. 2-22) 

(III. C. 2-23) 

(III. C. 2-24) 
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Figure 3-26 Pole-Zero Plot; Uncoupled Roll Channel Autopilot; 

Estimator Design; Discrete Closed Loop System; 
Circular Airframe 
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ROLL ANGLE, 
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Figure 3.27 Roll Angle vs Time; Uncoupled Roll Channel Autopilot; 

Estimator Design; Discrete Closed Loop System; 
Circular Airframe 
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Figure 3.28 Roll Angular vs Time; Uncoupled Roll Channel 

Autopilot; Estimator Design; Discrete Closed Loop 
System; Circular Airframe 
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Figure 3.29 Roll Tail Incidence vs Time; Uncoupled Roll Channel 

Autopilot; Estimator Design; Discrete Closed Loop 
System; Circular Airframe 
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Figure 3-30 Pole-Zero Plot; Uncoupled Roll Channel Autopilot; 

Simplified Estimator Design; Discrete Closed Loop 
System; Circular Airframe 
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Figure 3.31 Roll Angle vs Time; Uncoupled Roll Channel Autopilot; 

Simplified Estimator Design; Discrete Closed Loop 
System; Circular Airframe 
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Figure 3-32 Roll Angular Rate vs Time; Uncoupled Roll Channel 

Autopilot; Simplified Estimator Design; Discrete 
Closed Loop System; Circular Airframe 
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Figure 3-33 Roll Tail Incidence vs Time; Uncoupled Roll Channel 

Autopilot; Simplified Estimator Design; Discrete 
Closed Loop System; Circular Airframe 
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IV. MODERN CONTROL DESIGN AND ROBUSTNESS ANALYSIS OF 
COUPLED PITCH AND ROLL CHANNEL AUTOPILOT 
CIRCULAR AIRFRAME 

A. GENERAL 

The present chapter deals with the modern control design and 
robustness analysis of the discrete coupled pitch and roll 
channel autopilot for the circular airframe configuration. The 
continuous open loop coupled autopilot whose plant system and 
input matrices are presented in Appendix G is obtained by 
coupling the linear uncoupled pitch [Ref. 10] and roll (Table IV) 
channels. Then, utilizing analog- to-digital conversion by the 
aid of ORACLS program and for a sample period of 0.0125 seconds, 
the seventeenth-order discrete coupled system with matrices shown 
in Appendix H is formulated. Next, introducing the control-law 
and estimator designs were obtained and analyzed in terms of their 
transient responses and the application of the POPLAR design 
program [Ref. 7]. The POPLAR program is applied in order to 
employ singular value analysis and the use of an optimization 
routine to aid in pole placement control design of the above 
discussed linear multivariable systems. The robustness of the 
system is also considered by establishing singular value levels 
which correspond to multiloop gain and phase margins determined 
from the universal gain phase system diagram developed by Newsom 
and Mukhapadhyay at NASA Langley. 
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B. DISCRETE COUPLED STATE-FEEDBACK DESIGN 



1 . Design Approach and Analysis 

The discrete coupled state-f eedback designed autopilot is 
formulated by introducing into the original control system of 
Appendix H the following combined pitch [Ref. 10] and roll 
(III. B. 1-5) control-law gain vector. 



-c.oizs 


-QC279 


- 0.004-2 


a 029 


O. 004 L 


O 


0. 0001 


f./m 


’2.Q\62* 


0.0151 
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O 


0 


0 


O 


0 


O 




0 


0 


O 


0 


O 


0 


0 


0 


O 


0 


O 


O. 0004- 


0 


0. 0 1 79 


-0,01/9 


-1.499 


C. Go / 3 





(IV. B. 1-1 ) 

The pole-zero plot of Figure 4.1 indicates that the 
discrete coupled system is marginally stable, since the z-plane 



closed loop poles are: 

= 0. 1 3446+ jO . 0320974 
^ = 0. 13446-jO. 0320974 

^3 = 0.89699+jO. 0906457 
^ = 0.89699-jO. 0906457 
% = 0.955624+jO. 0299994 
= 0.953624-jO. 0299994 
^7 = 0.998202 
% = 1 

= 0.995608+j0. 0830713 
£,0 = 0.995608-jO. 0830713 



(IV. B. 1-2 ) 
(IV. B. 1-3) 
(IV. B. 1-4) 
(IV. B. 1-5) 
(IV. B. 1-6) 
(IV. B. 1-7) 
(IV. B. 1-8) 
(IV. B. 1-9) 
(IV. B. 1-10) 
(IV. B. 1-1 1 ) 
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Figure 4.1 Pole-Zero Plot; Coupled Pitch and Roll Channel 
Autopilot; State-Feedback Design; Discrete 
Closed Loop System; Circular Airframe 
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*2,, = 0.0993086 
£ a = 0.836555+jO. 306107 
2,3 = 0.836553-jO. 306107 
= 0.967788+jO. 0334452 
2,< = 0.967788-jO. 0334452 
£,6 = 0 . 89573+ J'O . 9381 33 
2 J? = 0.89573-jO. 938133 



(IV. B. 1-12) 
(IV. B. 1-13) 
(IV. B. 1-14) 
(IV. B. 1-15) 
(IV. B. 1-16) 
(IV. B. 1-17) 
(IV. B. 1-1 8) 



The time response plots of the roll angle, angular rate 
and tail incidence are presented in Figures 4.2 through 4.4. A 
close observation of these plots indicate that they are identical 
with those of the discrete uncoupled state-feedback design of the 
previous chapter. 

2 . Robustness Analysis 

Executing, the POPLAR design program of [Ref. 7] with 
inputs the data presented in Appendix I, the minimum additive 
input (MIN ADD IN SV) and output (SVADM0) singular values were 
computed from a frequency range from 0 to 200 rad/sec. 

Figures 4.5 and 4.6 which are plots of SVADM0 and MIN ADD 
IN SV versus frequency indicate that the discrete coupled state- 
feedback design is robust. It is noted that for very low 
frequencies the values of SVADM0 are above 0.8680. 

Finally, in terms of optimization results, the ordered 
computed eigenvalues are: 

Z, = 0.07717 (IV. B. 2-1 ) 

= 0.09855 (IV. B. 2-2) 

(IV. B. 2-3) 



Z 5 = 0.25923 
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Figure 4.2 Roll Angle vs Time; Coupled Pitch and Roll Channel 
Autopilot; State-Feedback Design; Discrete Closed 
Loop System; Circular Airframe 
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Figure 4.3 Roll Angular Rate vs Time; Coupled Pitch and Roll 
Channel Autopilot; State-Feedback Design; Discrete 
Closed Loop System; Circular Airframe 
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Figure 4.4 Roll Tail Incidence vs Time; Coupled Pitch and Roll 
Channel Autopilot; State-Feedback Design; Discrete 
Closed Loop System; Circular Airframe 
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Figure 4.5 SVADMO vs Frequency; Coupled Pitch and Roll 
Channel Autopilot; State-Feedback Design; 
Discrete Closed Loop System; Circular Airframe 
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Figure 4.6 MIN ADD IN SV vs Frequency; Coupled Pitch and Roll 
Channel Autopilot; State-Feedback Design; Discrete 
Closed Loop System; Circular Airframe 
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It 

«4 


0.069508 


(IV. B. 2-4) 


? _ - 


0.83707-jO. 3014 


(IV. B. 2-5 ) 


= 


0 . 83707+ jO .3014 


(IV. B. 2-6) 


it 

K. 


0.89844 


(IV. B. 2-7) 


= 


0.90064 


(IV. B. 2-8) 


It* 

11 


0.43996 


(IV. B. 2-9) 


2.o = 


0.46207 


(IV. B. 2-10) 


= 


0.96562-jO. 03401 


(IV. B. 2-1 1 ) 


£ u = 


0 . 96502+jO .03401 


(IV. B. 2-12) 


2,5 = 


0.99559-jO. 08307 


(IV. B. 2-13) 


2,4- = 


0.99559+jO. 08307 


(IV. B. 2-14) 


2,5 = 


0.99820 


(IV. B. 2-15) 


*6 = 


1 


(IV. B. 2-16) 


11 


1 .06546 


(IV. B. 2-17) 



C. DISCRETE COUPLED ESTIMATOR DESIGN 
1 . Design Approach and Analysis 

Following the same procedure as in the case of the state- 
feedback but for the coupled estimator gain vector the pole-zero 
and time response plots of the coupled estimator design are 
obtained. The time responses are again identical with those of 
the discrete uncoupled estimator design. 

2 . Robustness Analysis 

Executing the POPLAR design program of [Ref. 7], Figures 
4.7 and 4.8 are obtained which prove the robustness of the 
system. It is noted that for very low frequencies the values of 
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SVADMO are above 0.83976. Finally, in terms of optimization 



results the ordered computed eigenvalues are the same as in the 
coupled state- f eedback design. 
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Figure 4.7 SVADMO vs Frequency; Coupled Pitch and Roll 

Channel Autopilot; Estimator Design; Discrete 
Closed Loop System; Circular Airframe 
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Figure 4.8 MIN ADD IN SV vs Frequency; Coupled Pitch and Roll 
Channel Autopilot; Estimator Design; Discrete 
Closed Loop System; Circular Airframe 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The goal of the present thesis was the design and analysis of 
discrete lateral autopilots for application to BTT missiles. The 
following are the principal conclusions based on this work. 

1. The continuous and discrete classical designed autopilots 
were proved to have identical performances for the two 
lateral channels. 

2. The state-feedback and estimator autopilots were introduced 
as additional dynamic designs in order to implement control 
to the original system. Both designs, analyzed in terms of 
their transient responses, were found to meet the desired 
requirements . 

3. The simplified state-feedback and estimator designs reduced 
some of the returning gain loops, making the system 
simpler, without any significant effects on the system’s 
performance . 

4. The performance of the coupled pitch and roll channel 
autopilot was found to be satisfactory and the overall 
system proved to be robust. 



B. RECOMMENDATIONS 

In order to improve the simplicity of the overall system, 
more returning gain loops of the state-feedback and estimator 
design must be eliminated. A further investigation then must be 
conducted in order to examine if the performance of the resulting 
design remains unchanged. 
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APPENDIX A 



DESIGN REQUIREMENTS FOR UNCOUPLED AUTOPILOTS 

The requirements for the classical design method of the 
uncoupled channel autopilots [Ref. 6] are the following: 

1 . High Frequency Attentuation in Actuator Command Branch 

a. Uncoupled Yaw Channel 

It must be 15 db at 100 rad/sec and zero angle-of- 
attack and sideslip. This requirement will provide sufficient 
high frequency attenuation for 30 Hz actuator and for body 
bending modes when high frequency filters are added, but it 
limits the ability of the yaw autopilot to minimize sideslip 
angle . 

b. Uncoupled Roll Channel 

It must be > y 15 db at 100 rad/sec and zero angle-of- 
attack. This will provide sufficient high frequency attenuation 
for ^ 30 Hz actuator and for elastic modes when high frequency 
filters are added, but this requirement limits the speed of roll 
angle response. 

2 . Relative Stability for Both Lateral Channels 

Gain margins >6 db, phase margins ^ 30° with a goal of 
12 db and 50°. 

3 . Acceleration Time Response 

a. Uncoupled Yaw Channel 

(1) 63$ time constant of approximately 0.4 seconds. 
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(2) Overshoot ^ 10$. 

(3) Steady-state error need not be zero, 
b. Uncoupled Roll Channel 

(1) 63% time constant of 0.5 seconds. 

(2) Overshoot 10%. 

(3) Zero steady-state roll angle error. 



147 



APPENDIX B 



AERODYNAMIC DATA 



The overall classical design developed by Arrow [Ref. 6] was 
performed for the selected flight condition of Mach number 3-95 
at 60000 feet altitude. The corresponding aerodynamic data 
presented below were taken or derived from the ICAO standard 
atmosphere tables. 



TABLE VII 



AERODYNAMIC DATA ( M=3 . 95 , H=60kft) 



Temperature, T (°R) 



389.988 



Sonic Velocity, a (ft/sec) 



968.47 



Pressure, p (lb/ft^) 
Density, (lb-sec^/ft) 
Velocity, V (ft/sec) 



3825.4565 



0.0002238 



149.78 



Dynamic Pressure, q (lb/ft^ 



1637.145 
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APPENDIX C 



LINEARIZED AERODYNAMIC DERIVATIVES 

The linearized aerodynamic derivatives at the selected flight 
condition and about a zero trim angle-of-attack are provided 
below for both airframe configurations [Ref. 6]. 



TABLE VIII 

LINEARIZED AERODYNAMIC DERIVATIVES 
(M=3 . 95 , H=60kf t ; a=0°) 





CIRCULAR 


ELLIPTICAL 


C Y 


-0.065 


-0.043 




-0.025 


0.024 


C Y- 


0.021 


0.016 


Y 


-0.050 


-0.042 




0.031 


0.023 
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APPENDIX D 



MISSILE SIZING AND MASS PROPERTIES 

In order to provide a realistic missile based on the 
aerodynmamically tested configuration concepts [Ref. 5], the 
models were assumed to be 1/6-scale and the mass properties were 
developed corresponding to mass distribution which might be 
expected for missiles of this size. All the geometric and mass 
properties are presented in the following table. 

TABLE IX 

GEOMETRIC AND MASS PROPERTIES OF MISSILE CONFIGURATIONS 



Length, 1 (in) 

Max. Diameter (in) 

Max. Major Axis (in) 

Max. Minor Axis (in) 
c.g. distance from L.E. (in) 
Reference Length, d (ft) 
Reference Area, S (ft) 
Weight, W (lb) 

I xx (slug-ft 2 ) 

I zz (slug-ft 2 ) 



Circular 

168 

24 

100 . 8 ( 0.6 ) 
2 
TT 

2525 

40 

810 



Elliptical 

168 

41.57 

13.86 

100 . 8 ( 0.6 ) 

2 

TT 

2475 

40 

853 



150 



APPENDIX E 



PROGRAM LOGIC FOR APPLICATION OF ACKERMANN’S FORMULA 

The program logic for computing the control-law gain vector F 

T 

or the transpose of the estimator gain vector _K , taken from 
[Ref. 9 ], is given in the following table. 

TABLE 

PROGRAM LOGIC FOR APPLICATION OF ACKERMANN’S FORMULA 



1. Read in <1>. I'. T. and A,, the number of stales. 

2. Comment: first we will read in the desired pole locations in the \-plane. convert them 
to r-plane polynomial coefficients, and construct ai^). 

3. I *— identity matrix. A, x A. 

4. AI.IMIA — I 

5. A — 1 

6. If A - V. . go to step IS. 

7. Read in pole location k as </ - ;/». 

8. If h = 0. co to step 14. 

9. A, 2 e\pt</7 i cos hT 

10. A., *— expi 2(/7 i 

11. ALPII \ — ALPHA - (‘l» • <]> - * A.l) 

12. k — A - 2 

1 3. Go to step 6. 

14. .4 , — expo/ / ) 

15. ALPHA — ALPHA - <<l» - A, » I) 

16 . A — A - I 

17. Go to step 6. 

18. Comment: now we construct the controllability matrix. 

IV. C — I 

20. E - r 

21. A — I 

22. If A > A., go to step 28. 

23. Comment: replace column A of C bv E. 

24. C[ : k] - E 

25. A — A - I 

26. E — <|> • E 

27. Go to step 22. 

28. Comment: now solve for the control law. hist form d as the hist row of I. 

:</. E — I|.V, ; ] 

30. Solve PC = E for If. 

31 K = It ... ALPHA 

32. END 
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APPENDIX F 



ACKERMANN FORTRAN PROGRAM 



FILE: ACKERMAN „ATrIV Al 



IJC3 

INTEGER NS,K,IA, IER.I ,J,NI 

REAL Ph(20,20),GA(20,2J).AI(2C,20). ALPHA ( 20 , 20 ) , „K( AO ) 

REAL A2 ,t (20,20) ,0 ( 20, 23) ,3(20) , 6 MlJ,2 J) ,UM 20, 2 J) 

REAL A . 5 * T * R. A TEMP ( 20, 20 ) 

REAL «1 , T0120, 20) , T I C 23,20) , TS( 20,20) , 2(20,20), TM(20 ,20) 

C INPUT OF NS- NUMBER OF STATES 

C T= sample time 

PRINT , * ENTER NUMBER OF STATES.' 

REAU, NS 

PRINT, ’ ENTER SAMPLE TIME ' 

READ, T 

C INPUT JF P HI = A MATRIX 
00 100 1=1, NS 

OC 101 J=1,NS 
WRITE16.103) I,J 

103 FORMAT! AX, ' EN TEk PH( * , 12 , ' , ’ , 12, ' ) ' ) 

REAO.Ph ( I , J ) 

101 CONTINUE 
100 CONTINUE 
C INPUT OF G A=3 MATS IX 
00 200 1=1 , NS 

J = 1 

WRITE16.105) I,J 

105 FORMAT l AX , ' tNTEix GA ( * , 12 , ' , ' , 1 2 , ' ) ’ ) 

REAO.GA (1,0) 

200 CONTINUE 

C BUILDING THE IDENTITY MATRIX 
C Al PHA=C= I 

00 1 1=1, NS 
00 5 J = 1 » N S 
IF ( I ,E C.J ) GC TO 300 
A I ( I , J ) = J . 

ALPHA! I , J ) = 0 . 

Cll ,J) = 0. 

GC TU 5 

300 A It I , J ) =1 . 

ALPHA ( I ,J) = 1 . 

C ( I , J ) = 1. 

5 CONTINUE 

1 CONTINUE 
K = 1 

16 IF (K.GT.NS) GO TO IB 

C 

C INPUT OF DESIRED PULE:, LOCATION 
C 

WRITE (6, A3) 

A3 FORMAT!' ENTER REAL PART OF OESI RED POLE LOCATION.') 

READ, A 
WRITE! 6 »AA ) 

AA FORMAT!' ENTER IMAGINARY PART OF DESIRED POLE LOCATION.’) 

R EA 0 , o 

C BUILDING ALPHA (PH I) 

IF (8 . EJ.O ) GO TO IA 
A 1= -2.* EXP! **T> *COS<S*T ) 

A2= EXP ( 2.*A*T ) 

C EVALUATE P h I *P H I 

CALL VMULFF (PH, PH, NS, NS, NS, 20, 20, TM, 20, ICR) 

c 

OC 2 I = l, NS 

OU 6 J — 1 , NS 

TO(I , J ) =P H l I , J ) * A 1 
T I < I » J) =A I ( I, J) = A2 
TS( l, J) = TM( I,J)-t*TUCI»J)-t*TI ( I»J) 

6 CONTINUE 

2 CONTINUE 

CALL VMULFF ( ALP HA , T S , NS , NS , N S ,2 0 , 2 0 , A T EHP , 20 , I E R ) 

00 600 1=1, NS 
00 601 J= 1 , N S 

601 ALPHAl I ,J)=ATEMP(I ,J) 

6C0 CONTINUE 



ACKOOC 10 
AC ROOD 20 
ACK0C03D 
ACKOOOO 
ACKu0D50 
ACKCDObO 
AC K 00 0 TO 
AC KGuD 80 
ACK0C090 
ACKJ0100 
ACK0C1 10 
AC KOG 1 20 
ACKOul 30 
AC KQO 1 “tO 
AC K0U1 50 
AC KC C 1 60 
AC KOO 170 
AC KOO 1 80 
AC K Jo 1 90 
AC N002 00 
ACK002 10 
ACKJC220 
ACR0C230 
ACK002 AD 
ACK00250 
AC KC 02 60 
AC nC 02 70 
ACKJ0230 
ACK002 90 
Ai, KC J3 00 
AC K0D310 
ACK00320 
ACrDJ 33J 
ACKD03AD 
AC KO 03 50 
AC K0036D 
ACK003 70 
AC KC038J 
AC K00390 
AC KOoAOO 
ACK00A1D 
AC KC3A20 
ACK0CA30 
AC KODAAO 
AC KDUA 50 
AC KJOA60 
AC KO DA 70 
ACKGGABJ 
ACKO J A 90 
ACK0C5 JO 
AC KOO 5 1 0 
ACK00520 
AC K005 30 
AC KOO 5 AO 
AC KO 05 50 
ACK00560 
AC KOU 5 70 
ACK0U5 80 
ACK00590 
ACK00600 
ACK0U610 
AC K00620 
AC KDu630 
AC KOU6AO 
AC KJ0o50 
ACK00660 
AC K0U6 70 
AC K006 oO 
AC KO 06 90 
AC K007 00 
AC KO 07 10 
ACK0U720 
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FILE: ACKERMAN -iATFIV 



A l 







K= R+2 


a:\G07 30 






GO TO 16 


AC <0 u 7 ^ 0 




1 4 


A 1= E XP ( A*T) 


ACKC0750 






DO 3 1 = 1, NS 


A C n 0 G 7 6 0 






DO 7 J=l,NS 


AC/\0G770 






T 0 ( I , J) =Al l I , J) *A1 


AC ,\0 J7 30 




7 


T I ( I . J ) =PH( 1 , J J — TCI I ,J) 


AC <007 90 




3 


CCNT IhUt 


AC KG03 00 






CALL VMULFF l ALPHA, TI ,NS,NS.NS ,20,20 ,AT EMP.2 0, IcR) 


AC KOc 3 10 






DO 700 1=1, NS 


ACKOC 820 






00 701 J=1 , N S 


AC < 008 30 




701 


ALPHA ( I ,J)=ATtMPlI,J) 


AC<0o840 




700 


CON r INUE 


AC K 008 50 






K=K 1 


ACK00860 






GO TO 16 


ACK00670 




13 


00 9 1 = 1, NS 


AC K 008 80 






J = 1 


ACK00390 




9 


Et I , J) = GA( I , J 1 


AC K0C900 






K= 1 


AC <009 10 




22 


IF IK.GT.NS) GC TO 28 


ACK0092J 






DO 10 1=1, NS 


ACKOOy 30 






C ( I . K) = El i , 1 ) 


AC KO 09 40 




10 


CONT INU E 


aC K00950 






K=K *■ 1 


AC K 009 60 






CALL VMULFF IPH,E,NS,NS,1 . 2 0 , 20 , E , 20, I ER J 


AC K 00970 






GO TO 22 


AC K 0C9 80 




28 


00 11 J=1,NS 


AC<0j990 




11 


El 1 , J) =AI l NS.J ) 


ACkOIOCO 






R =0 . 


ACK010 10 






CALL L li I Nr l C ,2 0 »•', S , nS , R , 0 , 2 3 , S , «K. * I ER ) 


AC K 01020 


c 






AC. <01030 






CALL VMULFF IE,J,1 , NS ,NS , 20 . 20 , 8F , 20 , I ER ) 


ACKOiO^O 






CALL VMULFF IbF, ALPHA, 1 , NS , NS ,2 0 , 20 ,(J K, 20 , IE RJ 


AC <0 10 50 


c 






ACK0106C 






PRINT, * CGNTrUJL GAIN VECTOn. • 


AC <0 1 0 70 






WRITE ( 6,45) l UM 1, JJ ,J=1,NS ) 


A C <0 1 0 SO 




45 


FORMAT l 6F1 5 .4 ) 


A C < 0 1 0 9 0 






RETURN 


AC KO 1 1 00 






END 


ACKOlilO 


SENTRY 




AC K 01120 
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APPENDIX G: PLANT SYSTEM AND INPUT MATRICES OF CONTINUOUS COUPLED PITCH 

AND ROLL CHANNEL AUTOPILOT 
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APPENDIX I 



COUPLED STATE-FEEDBACK DESIGN INPUT DATA FOR POPLAR PROGRAM 
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